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We discuss and solve the transport without transit quantum paradox recently introduced in the 
context of the adiabatic transport of a single particle or a Bose-Einstein condensate between the 
two extreme traps of a triple-well potential. To this aim, we address the corresponding quantum 
dynamics in terms of Bohmian trajectories and show that transport always implies transit through 
the middle well, in full agreement with the quantum continuity equation. This adiabatic quantum 
transport presents a very counterintuitive effect: by slowing down the total time duration of the 
transport process, ultra-high Bohmian velocities are achieved such that, in the limit of perfect adi- 
abaticity, relativistic corrections are needed to properly address the transfer process while avoiding 
superluminal matter wave propagation. 
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Classical laws govern all of our daily physical situa- 
tions providing us with an intuition on what is physically 
feasible. It is not surprising then that a wide collection 
of quantum and relativistic phenomena escape from our 
'limited' classical intuition. Very often, the explanation 
of some puzzling non-classical effects has led to appar- 
ent contradictions that defy intuition such as, to cite 
only a few, the well known Schrodinger cat and Einstein- 
Podolsky-Rosen quantum paradoxes, or the tunnel and 
twin paradoxes in the relativistic realm. The detailed 
analysis of such paradoxes has allowed physicists to 
achieve a deeper understanding on the background of 
both quantum mechanics and relativity which, in turn, 
has strongly influenced modern research in quantum in- 
formation and cosmology, among other fields. 

A novel quantum paradox named transport without 
transit (TWT) plj was recently formulated in the context 
of the transport of a Bose-Einstein condensate (BEC) 
between the two extreme traps of a triple-well poten- 
tial by means of the matter wave analog of the quan- 
tum optical stimulated Raman adiabatic passage (STI- 
RAP) technique [2;. This spatial adiabatic transport 
technique has been studied for both single particles [3H5] 
and BECs [I] [7] and consists in adiabatically following an 
energy eigenstate of the system, the so-called spatial dark 
state, that only populates the vibrational ground states 
of the two outer wells and presents at all times a node in 
the central region. Thus, the matter wave is transferred 
between the outer traps, but maintaining an arbitrarily 
close to zero population in the middle region during the 
whole evolution [II H]. From this observation, M. Rab 
et al. have concluded that quantum-mechanically it is 
possible to transport matter directly from the left trap 
to the right one, without transiting the center region, 
and formulated the TWT paradox [T]: "Classically it is 
impossible to have transport without transit, i.e., if the 
points 1, 2, and 3 lie sequentially along a path then an 
object moving from 1 to 3 must, at some time, be lo- 
cated at 2. For a quantum particle in a three- well sys- 



tem it is possible to transport the particle between wells 1 
and 3 such that the probability of finding it at any time 
in the classical accessible state in well 2 is negligible." 
Clearly, quantum TWT is in contradiction with the con- 
tinuity equation that derives from both the single-particle 
Schrodinger equation and the Gross-Pitaevskii equation 
(GPE), that governs the ultracold BEC dynamics in the 
mean field approximation [5]. 

Bohmian mechanics [5] , while being equivalent to stan- 
dard quantum mechanics when averaged over the com- 
plete set of initial conditions, provides a very clear phys- 
ical picture of the quantum continuity equation. It is 
worth noting that average trajectories on a double slit ex- 
periment constructed from weak measurements of photon 
momenta have been experimentally observed recently, 
and are identical to the corresponding Bohmian trajecto- 
ries |10| . In this work we will numerically investigate the 
adiabatic transport of a BEC in a triple-well potential in 
terms of Bohmian trajectories to demonstrate that mat- 
ter wave transport always implies transit. We will show 
that by slowing down the transport process, Bohmian 
velocities in the middle well increase with no apparent 
limit. The appearance of superluminal velocities is "an 
artifact of using the nonrelativistic Schrodinger equa- 
tion" yjj, and thus in the hmit of perfect adiabaticity 
relativistic corrections would be needed to properly de- 
scribe the system dynamics. Ultimately, the origin of the 
TWT paradox is the incorrect use of the (nonrelativistic) 
Schrodinger equation in the adiabatic limit. 

The GPE governing the dynamics of a BEC reads: 
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where ^ is the wavefunction of the BEC of mass m, V 
is the trapping potential and g governs the non-linear 
interaction between the atoms [8]. For g = we re- 
cover the standard Schrodinger equation that deals with 
the dynamics of a single atom or a non- interacting BEC. 



In Bohmian mechanics [5], ip{f,t) guides an ensemble of 
Nt ^- oo trajectories, {n;[i]}, whose velocity is given by 

ri.[i] = ^;(f,i)|-^,^^[,], (2) 

\'ip(r,t)\'^ m 

S is the phase of the wavefunction and J is the probabil- 
ity current density. The trajectory ensemble reproduces 
the probability and current densities, 
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and thus, averages over the trajectory ensemble are 
equivalent to mean values obtained from the wavefunc- 
tion. Furthermore, by inserting the wavefunction in polar 
form, i.e., ijj = \ip\e'^^^'^, in the GPE and separating real 
and imaginary parts we are led to the equations: 
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Equation ([6]) is the so-called quantum Hamilton- Jacobi 
equation, analogous to its classical counterpart, but with 
an additional term, the quantum potential, that accounts 
for the quantum features of the system. Equation ([7]) 
is the quantum continuity equation for the condensate 
density that ensures that transport is always carried out 
in a continuous fashion and that Eqs. (|4]-l5]) are verified 
at all times [H] . 

It has been shown that the TWT paradox appears in 
both the ID and the full 3D cases [T]. Thus, to ease the 
computational burden we will work with the straight- 
forward reduction of Eqs. ^yj to one spatial dimension. 
Following Ref. [1] , we consider a triple- well potential con- 
sisting of a harmonic trap of frequency Wj, (and ground 
state width a~^ = y^h/mujz) where two Gaussian barri- 
ers of heights V12 and V23 have been added at positions 
=Fxo, giving a total potential of the form: 
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The barriers, created for instance with two blue-detuned 
lasers, will be considered initially large enough to inhibit 
tunneling such that three distinct traps are well defined. 
The eventual lowering of the barriers would allow the 
BEG, located initially in the ground state of the left trap, 
to tunnel to the other trapping regions. 

The adiabatic transport of the BEG from the left to 
the right trap is achieved by following an energy eigen- 
state of the system, the dark state, that only involves 




FIG. 1. (color online) Adiabatic BEG transport for g = Q and 
iOxtf — 5000. (a) Atomic probability density \'ip{x, t)\^ , barrier 
positions x — ±xo (dotted lines) and potential contour lines 
for V{x,t) = 10 hujx, 100 huj^ and 750 hw^ (from thinner to 
thicker), (b) Integrated populations Pi — L \ip{x,t)\'^dx with 
i — L, M, Rin the regions Cl = (—00, —xo), Cm = {—xo, xq) 
and Cr — (xo, 00). 



the ground states of the left and right traps, namely \L) 
and \R), and reads 1-0(6*)) — cos^jL) — sin6'|i?), where 
the mixing angle 9 is defined as tan0 — fli/fl2, being fli 
{^2) the tunneling rate between left and middle (mid- 
dle and right) traps [3]. By lowering the barriers in a 
counter- intuitive fashion, i.e., favoring first the tunneling 
between the middle and right traps and then the tunnel- 
ing between the left and middle traps, \D{9)) transforms 
from \L) to \R) [3\. With this in mind, the temporal 
profile of the potential barriers height is taken to be: 
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being t^ the time delay between the pulsing of V23{t) 
and Vi2{t), giving a total adiabatic transport time T = 
tp+td- For the numerical simulations, we take parameter 
values similar to Ref. [L: T/min = ^fy^x, I^nax — ^Q^fy^x, 
a — 0.16q:"^ and xq = 0.48a~^. We will simulate the 
process for different values of g to study the effect of the 
nonlinear interaction. We will also analyze the effects of 
slowing down the process, by varying the pulsing time 
ip, even though for tOxtp < 1000 the dynamics is not 
adiabatic enough and the transport breaks down. The 
delay id is in all cases taken to be 0.15tp. 




FIG. 2. (color online) (a) Bohmian trajectories Xk[t\ as- 
sociated to the evolution shown in Fig. Ilja) (black lines), 
atomic mean position (solid line), dark state node position 
a;"(f) (dashed line), and barrier positions x — ±xo (dotted 
lines), (b) Bohmian velocities Xk[t] of the trajectories in (a) 
(black lines), velocity of the atomic mean position (solid line), 
and Bohmian velocity at the node position (dashed line). 



We show in Fig. flla) the probability distribution dur- 
ing the BEC transport with the barrier heights from 
Eqs. (|9- 10). Indeed, when the potential barriers are low- 
ered it appears as if the condensate is vanishing from the 
left region and appearing in the right region without tran- 
siting the space between the barriers. This is validated 
by Fig. Il|b), that shows that while the BEC population 
is completely transfered from the left to the right region, 
the middle region remains unpopulated the entire time, 
due to the fact that the dark state that the system is 
(ideally) following presents a node in this region. 

To elucidate the transport process we have computed 
the Bohmian trajectories associated to this evolution, 
that are shown in Fig. l2ja). The trajectories, initially 
distributed according to\'il){x,t — 0)p, follow the wave- 
function at all times, starting in the left trap and end- 
ing in the right one. All the trajectories cross the mid- 
dle region, demonstrating the non-existence of transport 
without transit, but in order to keep a low population 
of the middle region by minimizing the time spent there, 
they are forced to increase their velocity [6]. We com- 
pute the dark state node position, a;"(i), as the position 
in the central region where the population is minimum, 
and see that the trajectories accelerate when they get 
closer to it, reaching velocities much larger than the mean 
wavepacket velocity, see Fig. WCb). Note that each tra- 
jectory velocity peaks at a different time, t^, depending 



on its initial position inside the wavefunction, but all 
trajectories achieve similar maximum velocities, t^ cor- 
responds to the time at which the fc-th trajectory crosses 
the dark state node, i.e., Xkif^] = 2:"(^"), since the trajec- 
tories peak velocities correspond to the Bohmian velocity 
at the node, see dashed line in Fig. [2Fb). 

Since the transport is performed in a finite time, the 
dynamic state ip{x,t) will not be exactly equal to the dark 
state and will present some small population from other 
states giving a non-zero atomic population at the node 
position during the evolution. If the evolution is made 
more adiabatic by performing the process more slowly, 
the dynamic state will remain closer to the dark state 
and the node population will be smaller. This will result 
in an increase of the trajectories velocity at the node (cf. 
Eq. (pi)), presenting a surprising effect: by slowing down 
the total matter wave STIRAP sequence it is possible 
to achieve sudden trajectory accelerations yielding ultra- 
high velocities. 

In order to obtain a quantitative argument, we com- 
pute the ensemble average of the velocity that the trajec- 
tories achieve at the node, {v"^'^^). By using Eqs. pHsl) 
and taking into account that Xk[t] > and i"(t) < (cf. 
Fig. I2|a)) such average can be rewritten as: 
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The second term in the last integral of Eq. ( |Tl| ) is orders 
of magnitude smaller than the first one (and decreases 
as we increase the adiabaticity) since i" is much smaller 
than the peak velocities. 

We plot in Fig. [3] the result of various simulations for 
different times ip and different interaction strengths g. 
Variations in g slightly change the values of the calcu- 
lated quantities, but the overall behavior remains the 
same. Thus, for clarity we have performed fits for the 
.9 = data only. Figure Isla) shows that the mean ve- 
locity that the trajectories reach when passing through 
the node grows linearly with tp. Note that all trajec- 
tories accelerate abruptly, but at different times, and 
then, the maximum value that the mean atomic veloc- 
ity, max [d(a;)/cft], takes during the evolution decreases 
inversely proportional to tp, see Fig. pTb), remaining or- 
ders of magnitude below {v™'^^). There is no apparent 
limit to the trajectory velocities in the middle region as 
we approach the limit of perfect adiabaticity, and at some 
point they might surpass the speed of light. Bohmian 
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FIG. 3. (color online) Results for different simulations with 
different values of ip and g. (a) Mean trajectory velocity at 
the node, (b) Maximum mean atomic velocity, (c) Maximum 
probability density at the node. Solid lines are, respectively, 
a linear fit, a tp ^ fit and a tp^ fit for the g = data. 



velocities calculated from the Dirac equation are always 
bounded by the speed of light ^llj, implying that rela- 
tivistic corrections to the Schrodinger equation would be 
needed to properly address the transport while avoiding 
superluminal matter wave propagation in the adiabatic 
limit. It is surprising that the Schrodinger equation (or 
the GPE) cease to be valid and one should consider such 
corrections in the limit where the process is performed 
"infinitely" slow. It is also striking that, at variance with 
the usual tunneling problems where such high-velocities 
occur at the energetically-forbidden regions, here they 
appear in a minimum of potential. 



Figure Islc) shows that the maximum population of 
the node during the evolution decreases inversely pro- 
portional to tp. As we already mentioned, this non-zero 
population appears due to an imperfect following of the 
dark state, and its decrease is what makes the high ve- 
locities in the middle region to grow, as a classical fluid 
is forced to accelerate when traveling through a narrow- 
ing conduct. Note that a perfect node would forbid the 
transport, and thus a non-zero population is needed for 
the transport to take place. We have computed the in- 
tegrated flux through the node, /^ J{x'^{t),t)dt, and as 
expected, found it to be 1 in all cases, meaning that the 
entire wavefunction that is transported from the left to 
the right region is transiting through the (quasi)node in 
the middle region. 

In summary, by using the Bohmian formalism of quan- 
tum mechanics applied to the adiabatic transport of a 
BEG in a triple well we have shown that quantum trans- 
port always implies transit. The absence of atomic pop- 
ulation in the middle region is explained by an increase 
of the trajectories velocities in the vicinity of the node of 
the followed dark state. As the process is performed more 
adiabatically, the trajectories velocities increase with no 
apparent limit, and thus relativistic corrections should be 
taken into account to describe correctly the dynamics of 
the system. It would be interesting to investigate if the 
Doppler shift in the light absorption due to the atomic 
velocities could be used to (weakly) measure the atomic 
momenta and study some of the features described in this 
paper in similar lines as in Ref. [10,. Moreover, during the 
adiabatic transport, the velocity change is very abrupt, 
leading to very high accelerations (and decelerations) . It 
remains an open problem to investigate whether the mat- 
ter wave STIRAP technique for a charged particle could 
lead to the emission of radiation. 
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